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Abstract 

The reflection of a three-dimensional vectorial Maxwell- Gaussian beam by a planar surface is 
studied. The surface is characterized by its complex reflection coefficients rs(k) and rp(k) for 
TE and TM electromagnetic plane waves of wavevector k, respectively. The field impinging upon 
the reflecting surface is modeled as a quasi-monochromatic fundamental Gaussian beam suitably 
modified in order to satisfy Maxwell equations (Maxwell-Gaussian beam). Analytical expressions, 
correct up to the second order in a perturbation expansion, are given for the refiected electric and 
magnetic field, respectively. We found that first order terms in the perturbation expansion account 
for a longitudinal shift (Goos-Hanchen effect) of the whole reflected beam, while second order terms 
modifies the transverse shape of the beam which is, at this order, no longer cylindrically symmetric. 



I. INTRODUCTION 



Plane waves of the form wexp(ik ■ r — iut), (with u ■ k = and uj = c|k|), are solutions 
of the Maxwell equations that show a serious problem: They are physically impossible since 
posses an infinite amount of energy. However, they are very easy to handle and, thus, widely 
used in the physics community. Moreover, an electromagnetic field of a given arbitrary 
shape, can always be written as a linear superposition of plane waves, hence, for example, 
it is enough to know how a plane wave propagate across some medium, to know how the 
whole field does. For this reasons, even in advanced textbooks^, the reflection coefficients 
associated to the interface between two media, are calculated in terms of the amplitudes of 
incident and reflected plane waves. However, since in our real world plane waves of infinite 
transverse extension do not exist, but only finite-transverse-size light beams, non-specular 
reflection effects are expected^'^'^'^'^i^ and, actually, occur. The most known effects are 
the Goos-Hanchen^i^ii^iii and the Imbert-Fedorov^^ii^ longitudinal and transverse shifts^, 
respectively. 

In these Notes we study in a detailed and didactic manner the non-specular effects oc- 
curring when a Maxwell-Gaussian beam^^ impinges upon an arbitrary planar surface char- 
acterized by its complex reflection coefficients^. 

II. GEOMETRIC REFLECTION 

Before discussing the physical process of reflection of light by a planar surface, let us 
consider some general characteristic of geometric reflection by an ideal planar mirror (for 
a short and simple introduction to mirror symmetry applied to electromagnetism see, e.g., 
ref.— ). Let K = (Oxyz) be a Cartesian reference frame whose axes x,y,z are specified 
by the three unit basis vectors {sx, By, e^}, respectively, and let 2; = be the equation of 
the reflecting planar surface. The geometry of the problem at hand is illustrated in fig. 1. 
We define geometric reflection (or mirror symmetry with respect to the plane z = 0) in an 
operative fashion, as follows. If a point P has coordinates r = e^x + By-y + e^z in K, then 
the mirror image point P has coordinates r = e^x + Cyy — e^z. Thus, the mirror image of 
a scalar field (j){x, y, z), is simply a new field (j){x, y, z) defined as 
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FIG. 1: Illustrating the geometric reflection (or mirror symmetry) of a right-handed circularly 
(RHC) polarized plane wave by an ideal mirror. The grey plane of equation z = is the ideal 
mirror surface. 

Next, we must consider geometric reflection of vector fields. To this end we study, without 
loss of generality, the reflection of a circularly polarized plane wave with wave vector k, 
impinging upon the xy-plane from 2; < 0, as shown in Fig. 1. Geometric reflection changes 
k in k and transforms the image of a right hand in the image of a left hand. Thus, the 
mirror image of a right-hand circularly polarized plane wave, must be a left-hand circularly 
polarized plane wave, and vice versa. These statements can be straightforwardly transformed 
in formulas as follows: Let us write the incident plane wave of unit amplitude as 

where cc;(k) = c|k| and the orthogonal basis vectors {ei(k)}?^^ are defined as in Fig. 1 and 
Appendix A. According to the definition of geometric refiection given above, the refiected 
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field A (r, t) is still a plane wave of unit amplitude with equation 

A-^(r, t) = ^i(k)-^^^2(k) ^ ~ _ ^ _ ^^^^^^^ 

where k = e^k^ + e.yky — e^k^ lies on the plane of incidence defined as the common plane of 
k = e^kx + Gyky + Bzkz and e^. Note that k ■ k = k ■ k, therefore |k| = |k| ^ ^^(k) = a;(k). 

It is easy to generalize the results above to the case of an arbitrary incident field. By 
definition, such a field can always be written as a plane waves expansion of the form: 

A^°^(r,t) = j [ei(k)ai(k) + e2(k)a2(k)] exp [ik ■ r - icj(k)t] A^k. (4) 

Since from Eqs. (EI)-® we know how each plane wave making A™'^(r,t) transforms under 
geometric reflection, we can write at once 

A^^f(r,t) = y"[ei(k)ai(k)-e2(k)a2(k)]exp[ik-r-ia;(k)t]d3fc 

= ^2 ex{k) rxax{k) exp\ik -r - iuj(k)t]d^k, (5) 

A=l ^ 

where we have defined the reflection coefficients 

rx= < (6) 
[ -1, A = 2, 

and we have used the property k ■ r = k ■ r . Note that according to Eqs. ©-Q, the 
amplitudes ax(k) do not change by geometric reflection, that is aA(k)-\> aA(k). 



III. PHYSICAL REFLECTION 

Now, let us consider the case of reflection by an actual physical planar surface charac- 
terized by the reflection amplitude coefficients (with respect to the vacuum) r<j(k), rp(k) for 
s-polarized waves (or TE, that is plane waves with the electric field orthogonal to the plane 
of incidence), and p-polarized waves (or TM, that is plane waves with the electric field par- 
allel to the plane of incidence), respectively. From our choice (see appendix A) for the basis 
vectors {ei{k)}f^i, with ei(k) and 62 (k) parallel and orthogonal with respect to the plane 
of incidence, respectively, it follows that rp(k) = ri(k), rgik) = r2(k). For an homogeneous 
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medium with complex-valued dielectric constant e and refractive index^ h = \fk = n + in, 



(n, K G M) , we have 



where k = e^kx + Byky + Bzkz, and = exkx + eyky + Czkmz is the wave vector of the plane 
wave transmitted into the medium, as given by the Snell law, and 



k^z = ^ikl + {e-l){kl + kl), (9) 

is the (generally complex- valued) 2;-component of k^. 

Once we know the reflection coefficients {'"A(k)}|=i associated to a single plane wave with 
wave vector k, we can easily determine the behavior under reflection of an arbitrary field 
just by letting r\ — > rA(k) in Eq. ([5]): 

A'^\r,t) = J2 eA(k)rA(k)aA(k)exp[ik-r-ia;(k)t]d3A;, (10) 

A=l 

Equation ( fTOl) is perfectly general and, therefore, of limited usefulness. However, much 
additional work can be done if we consider the actual experimental situation where the 
incident field is a quasi-monochromatic narrow beam directed along kg with central frequency 
Uq = c|ko|. Such a beam can be represented by an envelope vector field ^''"'^(r, t) modulating 
a carrier plane wave with wave vector ko and frequency ujq: 

Ai-(r, t) = *"-(r, t) exp [z(ko • r - cuot)] • (11) 

The envelope ^'™'^(r,it) is easily determined by rewriting Eq. (fTT!) as 

A^°^(r,t) = exp [z(ko ■ r - cuot)] 

2 „ 

/ |eA(k)aA(k)exp [z(k-ko)-r] 

A=l ^ 

X exp [-t{uj{k) - iUo)t] jd^fc. (12) 



^ Warning: Note that the present definition of n is different from the one given in Chap. Xfff of the book 
by Born and Wolf, where the authors write n = n(l + zk) 
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The very same procedure can be executed for the reflected beam obtaining an envelope 
vector field ^^^^{r,t) modulating a carrier plane wave with wave vector ko and frequency 



A^^f (r, t) = *'""^(r, t) exp z(ko • r - c^o^) 



where 



*-^(r,t) = J2 {eA(k)r;,(k)a;,(k)exp[z(k-ko)-r 

A=l 



exp [-i[uj(k) -uJo)t] \d^k 



Since (k — ko)-r =(k — ko)-r we can write 

*-^(r,t) = J2 {eA(k)r,(k)a;,(k)exp[2(k-ko)-?] 

A=l ^ 

X exp [-^(^^(k) - coo)t] jd^/c 
2 „ _ 

= ^ |eA(ko + q)rA(ko + q)aA(ko + q)exp(iq-r) 

A=l 



X exp [—i5uj(ko,q)t] \d^q 



(13) 



(14) 



(15) 



were in the second line we have changed the variables of integration according to: 



k -> q = k - ko, d^k = d^g, 



(16) 



and we have defined 



5uj{ko, q) = u;(ko + q) - ujq. 



(17) 



By hypothesis, we are considering a narrow beam, therefore we expect that the difference 
vector q = k — ko would be small as compared to ko: |q|/fco -C 1^^. Then, we can Taylor 
expand the reflection coefficients rA(ko + q) around q = 0. To this end, first we rewrite 
rA(ko + q) as 

rA(ko + q) = |rA(ko + q)|exp[i(/)A(ko + q)] 

= i?A(ko + q) exp [iM^o + q)] , (18) 
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FIG. 2: Illustrating the geometric meaning of the change of integration variables k ^ q = k — kg. 

and then we make a Taylor expansion for both the amplitudes -RA(ko + q) and the phases 
0>,(ko + q), separately. For the latter we have, up to second order terms^^: 



^(ko + q) 



.(ko) + 



dki 



1 5Va 



k=ko 2 dkidkj 



QiQj + ■■■ 



k=ko 



= (pxC^o) + y^A ■ q + -q • ^Aq + • • • 



(19) 



where summation over repeated indices is understood. The "displacement" vectors ip\ 
cpA(ko) and the symmetric 3x3 tensors J^x = jFx(ko) are defined as 



<^Uko) 



•n^(ko) 



dkidkj 



k=ko 



.(k) 



(20) 
(21) 



k=ko 



G {1,2,3}), respectively. A similar Taylor expansion can be also performed for the 
amplitudes -RA(ko + q) after their exponentiation: 



Rxiko + q) = exp [lni?A(ko + q)] 



(22) 
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where 



ln/?A(ko + q) = \nRx(ko) + qi 



1 d^R 



1 dRx 
Rxik) dh 

A 



k=ko 

1 dRx dRx 



Rx(k)dkidk-j Rl{k) dh dkj 



k=ko 



1 



= lnRx(ko) + px ■ q + -q- TZxq 



(23) 
(24) 



The "amphtude modulation" vectors px = PxiMo) and the symmetric 3x3 tensors Tlx 
7^A(ko) are defined as 

1 



P\(ko) 
^r(ko) 



i?A(ko) 
1 



k=kn 



i?A(k, 



Oj 



dkidkj 



Rx{k) 



-pUko)pi(ko 



(25) 
(26) 



k=k| 







respectively. It is important to stress that all partial derivatives in Eqs. fll9p - (|2S]) are 
understood to be whole-partial derivatives^^, that is, if / = /(k, Ci;(k)) is a smooth function 
that depends both explicitly upon k via /(k, u;(k)), and implicitly upon k via u;(k), then 
with d/dki we mean 



|-/(k..(k))« 



d duj d 
dki dki duj 



/(k,a;(k)). 



(27) 



Now, by combining Eq. (fT8|) with Eq. (fT9l) and Eq. (l23l) . we obtain 



rA(ko + q) = Rx(}^o) exp ( Pa ■ q + -q • 7^Aq + 



X exp ( #A(ko) + ^V^A ■ q + -q ■ J^x^. + • • • 



rA(ko) exp 



(pa + iipx) ■ q + ^q ■ ("^A + «^A)q + • • • 



(28) 



where rA(ko) = i?A(ko) exp[i0A(ko)]. Finally, we substitute this expression into Eq. (fT5|) to 
obtain 

W^^{r,t) = J2 {eA(ko + q)rA(ko + q)aA(ko + q)exp(iq-r) 

A=l 

X exp [—iSuj{ko,q)t] \d^q 



^rA(ko) / eA(ko + q)aA(ko + q)exp[-i5t^(ko,q)t] 



A=l 

X exp < 



(r + c^A - ^Pa) + ^(J^a - i'R-x)q + ■ 



d^q. (29) 



This is our final expression. One can see that the whole beam undergoes both a linear 
— ip\ and non-linear {T\ — iR,\)(\ complex displacement. 



A. Reflection of vector Gaussian beams 

In this subsection, we consider an incident vector field tailored as a Gaussian beam. From 
Appendix D, we know that such a field can be described as 

2 



E'^^^lr.t) = J^i2j e,ik)a,{k)exp[z{k.r-ut)]d'k, 



(30) 



where 



a,{k) = a{k)e,{k)-f^a{k)f,{k). 



(31) 



The complex-valued unit vector / specifies the orientation of the polarizer that selects the 
polarization of the incident beam. Equation ( l30l) furnishes an expression valid for the electric 
field only. However, it is not difficult to see that we can also write the magnetic field as 

2 



B'-(r, t) = g I e,{k)b,{k) exp[^(k ■ r - ut)]d\ 



were 



J2ex{k)b,{k) = - 



A=l 



k X ^eA(k)aA(k) 



A=l 



(32) 



(33) 



Explicitly, we have 



h{k) = — a(k)/2(k), 

c 

b2{k) = -a(k)/i(k). 



(34) 
(35) 



were c is the speed of light in vacuum. Thus, if with A™'^(r, t) we denotes either the electric 
or the magnetic field, we can write it as 

(36) 



^'"^"'^^ = J^iZ J eMg,{k)exp[t{k.r-ut)]d'k 

LI — 1 



with 



9f.ik) = a(k) ^£^^e^(k) ■ /, 



(37) 



where we have defined the 2x2 matrix £ = [C,fj,v] as 



















= < 


















for the electric field, 
for the magnetic field. 



(38) 



Since we want to express A"^'^(r, t) in terms of the Fourier transform F(p) of the Gaussian 
beam /(r), we follows the procedure explained in Appendix C and make the change of 
variables k ^ q, where k = ko + q, thus obtaining: 

A'"'(r,t) = (2^)3/2 / e^(ko + q)^/.(ko + q)exp[i(ko + q) •r-iu;(ko + q)t]d3g 



^^(ko-r-tJot) 



(2vr)3/^ ^ 



^ / e^(ko + q)5f^(ko + q)exp[iq-r-i5u;(ko,q)t]d^g, (39) 



where 5ti;(ko, q) is defined by Eq. flTTl) . Under refiection this equation changes according 
the following rules: r — r, for the scalar part of the field, and 

2 2 



^ e^(ko + q)£^, [e,(k) ■ /] ^ X] ^^('^o + ^^^^^ ^^^^^^ ' -^l '"^^'^o + q) 
for the vector part of the field. Thus, from Eq. fl39l) it readily follows 



(40) 



A-'(r,t) 



^ / e^(ko + q)r^(ko + q)5f^(ko + q) 



(2vr)3/^ ^ 
X exp[iq ■ r — i6uj{kQ, q)t]d'^g 



(ko + q)r^(ko + q) 



(27r)3/2 

xa(ko + q) exp[zq ■ r — ^^^^(ko, q)t]d^g ^ ■ / 
= ^'■^f(r,t)-/ 

where Ai^'^^{r,t) is a 3 x 3 space-time dependent matrix. From Eq. flC27p it follows that 

a(ko + q) = F{Dq) exp[-i(ko + q) ■ fq]. 



(41) 
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which can be used in Eq. ( BTj) to obtain 

gj[ko-(r-ro)-<^0<] f ~ ^ 

-M"^(r, t) = j^^^ 5J l^, J e^(ko + q)e,(ko + q)r,(ko + q) 

xF{Dq) exp[zq ■ (r — Fq) — ^^^^(ko, q)t]d^g 



^i{koz"-uot) ^ 



(27r)3/2 



e^(ko + q)e^(ko + q)r^(ko + q) 



xF{Dq) exp[iL>q • Xx" - i5uj(ko, q)t]d^q, 



(42) 



where we have used Eqs. ( ]C49tiC50lC53P to rewrite ko ■ (r — tq) = ko ■ (r — fq) = koz", and 



q - (r -To) = Dq - D{r - Fq) = Dq ■ XD{r - ?o) = ^q ■ ^x". 



(43) 

At this point, Eq. is exact but not very usefuL However, we can exploit the hypothesis 
that the impinging hght beam is modeled as a narrow Gaussian beam, and make a Taylor 
expansion about k = ko (or, equivalently, q = 0), of all relevant quantities in Eq. fj42|) . In 
practice, we substitute into Eq. fH2|) . the expansion (128|) : 

q - {pu + i^u) + ^q • {'R-u + iJ^u) q 



r^(ko + q) = r^(ko)exp 



= r^(ko) exp ( ^ • + ^q ■ W^q 



r^(ko)exp 



■I] 



(44) 



where we have chosen to truncate the Taylor series of the amplitudes -Rjy(ko + q) at the 
second order terms, and we have defined v^, = A;o(p!/ + i^u), and = k^ilZy + iJ-'y). 
In practice, we have kept all the second order terms with the aim to deal with Gaussian 
integrals which are analytically integrable. Moreover, we use the results of Appendix F to 
write the second-order Taylor expansions of the diadic forms e^(ko + q)ei,(ko + q) as 



,(ko + q)e,(ko + q) ^ E^, + |^E« + i^M^^^^ ■ 



(45) 



By substituting Eqs. ( 14411451) into Eq. (142]), we obtain 

i{koz"-ujot) ^ 



(27r)3/2 



fJ,,U=l 



E + _E^*^ -\ 



X exp ( ^vi + ^qiqjWJ ) F{Dq) exp[tDq ■ Xx" - t6tuiko, q)t]d3g. (46) 
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Since, from Eq. (ICllI) we know that 



exp 



2h 



5 P3 + 



vl+vl 

2ko 



(47) 



it is convenient to perform a change of integration variables into Eq. (I46l) by putting q 
p = Dq ^ q = D^p, thus obtaining 



J{koz"~uJot) 

(27r)3/2 



kn 



2^0 



X exp 



PaDar+P^^D:^, ^^^ ^ J_(p^/)^. _^ p-iD-ii){p pD pj + P3D3j)WJ 
Kn ZKn 



X 



F(p) exp[i(p„X^ja;J + pg^si^) ~ *'^^(ko, £'^p)t]d3p, 



(48) 



where we have displayed the variable p^ in order to integrate with respect to it by exploiting 
the Dirac delta present in the expression of F{p). Such integration yields the substitution 

1 



Pa 



2ko 



PaPfsS, 



(49) 



everywhere in Eq. (HHl) . In addition, from Eq. (I49l) and Eq. (1C33I) . it immediately follows: 

^ik,,D-p)^^(dp^]\o, (50) 



where the last approximate equality follows from our choice to keep expansion terms up 
to the second order in the phase terms. This result is important because it tells us that 
up to the second order terms, the beam is monochromatic. Therefore, we are enabled 
to use, for the calculation of the Poynting vector S(r), of the cycle-average theorem to 
obtain S(r) oc Re [E(r,t) x B*(r,t)]; we shall use soon this result. Within the same level of 
approximation, we can thus rewrite Eq. fjlSj) after integration with respect to ps, as: 

i{koz" -uot) 2 /■ p ^ 



\Pa_Pll 

2 ka ko 



(27r)V^L 
kQ 



X exp I — [DaiVl + ikoXaiX ' 

PaPi3 [Da,DpjUi^ - 6a(3 {D^ivi + ik^z" + k^L)] \ d^p, (51) 



2k Q 
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where we have used Eq. ( ]C45P to write X^ix"^ = S^i^—lY^^x"^ = x 



~^ z" . If we introduce 
the dimensionless variables {Pq, = Pq//^o}q=1) then we can rewrite Eq. (|5T|) as 



ZTT ^ — ^ 



E 



1 



2 

X exp 



■-Pa(P;/ )al3Pf3 + Pab 



d^P, 



(52) 



where we have defined the two two-dimensional vectors h,^ {u = 1,2), and the two symmetric 
2x2 matrices B^, whose elements are 



= iko[D{ip^-ip^)+Xx."] 



5ap (&.3 + hL) - ikl[D{T, - in,)D'^]^ 



(53) 



and a, {3 G {1,2}. Note that in the case of geometric refiection, the expressions above reduces 
to = ikoXx.", and {B^^)a/3 = ik^dap {z" — iL), respectively. Moreover, it is important to 
note that the off-diagonal terms of the matrix B, which change the profile of the beam, arise 
because of the second order terms ^ Tv — ilZu] while fist order terms — ipy simply 

produce a shift. Equation (1521) can be analytically integrated by using the following three 
basic formulas for Gaussian integrals^i; 

2 2 



G,. 



exp 



a=l 



d^P 



27r(det B^fl'^ exp 



I 2 



a,f3=l 



(54) 



r 2 2 

f P.exp PaiB;')a(3Pp + J2P»^ 



d^P 



dG, 
dhun 



Gy {bi,a{Bp 



Gu {Byhy)^ . 



(55) 



G,, 



P^Pr exp 



2 2 
\ Y P.{B-')c.pPp + YPc.K 



Q!,/3 = l 



a=l 



d'P 



d'G, 



Gy [(P;/)ar + {Ka{.By) a^){hyp{By) Pt-)] 



dbadbr 
Gy [By + P,(b,b,)P, 



(56) 
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By using the equations above into Eq. (152|) . we easily obtain our final result: 

k r ^ 

fl,U = l 

h T ^ 

[5, + B,{KK)B,]^p (5«/3/^3.E« - Z},,D^,E(:^)) }, (57) 

with z,j G {1,2,3}. From this equation we can calculate analytically both the electric and 
magnetic fields and, consequently, the Poynting vector. From the latter we have automati- 
cally the energy fiux which is the quantity experimentally accessible. 

B. First order Taylor expansion 



Qip Qij) Qip 

ipikox + Qx, koy + Qy, koz + Qz) = ^{kox, koy, koz) + { ^Ix-q^ + ^ly^ + 9^ 



dk,, 



ifikox, koy, koz) + q • VkV'|k=ko + ' • • 



+ 



k=ko 



(58) 



where 



_ d d d 

Vk = eiT^ + 62771- + 63 



' dkn, 



Id Id d 

^ kdd ksm9 d(f) ^ dk 
/I N 1 /, \ 1 d . d 



If il){kxi ky, kz) = ip{0) then 

^'(ko + q) = V^(ko) 



1 dip 



ei(ko) ■ q + . 



e=eo 



ip{ko) + — [q^ cos e'o - Qz sm fc^o) ^ 



+ . . . 



(59) 



(60) 



14 



27r 



f2 



/iGi + y" cot ^o^r^ (Gi + PG2) 
z" - iL 



z" - iL 



hpG^ - y" cot ^0^7^ (^1 + 



27rc 



e*('=°^"-'^°*Vi(ko)- 



- 60 



f2pG2 - y" cot ^o^r^ (Gi + pGa) 

hC^ + y"cot9o^/^ (Gi + pGs) 

f2pG2{x"-v',)-y"hG, 
z" - iL 



APPENDIX A: COMPLETE BASES IN W 

Let K = iOxyz) be a Cartesian reference frame, and let k denotes a vector pointing along 
the direction fc : k = A;fc, where k = |k| and |fc| = 1. Moreover, let {e^, Cy, e^} be three 
mutually orthogonal unit vectors pointing along the axes x, y and z. respectively. Given 
k and Cz, we can built a complete basis {ei(k), e2(k), e3(k)} in by using the following 
recipe: 

{cz X k) X k 



ei(k) 



e2(k) = 
e3(k) = 



\{ez X k) X k| ' 
X k 



|e^ X k| ' 



Ikl' 



where " x " denotes the ordinary vector product in 



(a x h)i = eijkajbk, {i,j,k e {1,2,3}), 



(Al) 
(A2) 
(A3) 

(A4) 



and Cijk is the completely antisymmetric Levi-Civita tensor, such that eijk — +1 or —1 
according to whether the indices i,j,k are an even or odd permutation of the ordered set 
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{1, 2, 3}, and e^-fc = when at least two indices are equal. By definition, the three real- valued 
unit vectors {ej(k)}?^]^ form an orthogonal 

e,(k)-e,(k) = 5,^-, (A5) 

3 

ei(k) X ej(k) = ^eijfeefc(k), (A6) 

k=l 

G {1,2,3}) and complete basis in M^: 

3 

5^e,(k)e,(k) = /, (AT) 

i=l 

where / denotes the 3x3 identity matrix, and we have used the symbol ab to denote 
the diadic product of the two vectors a and b, respectively, that generates a matrix whose 
elements are: (ab)jj = aibj, = 1,2,3). Since e3(k)e3(k) = kk/k'^, from Eq. (1A7P it 
readily follows that 

^e,(k)e,(k)=/-^ (A8) 

A=l 

= I-C(k) (A9) 

^r(k), (Aio) 

where with C{k) and T(k) we denoted the longitudinal and the transverse projectors, re- 
spectively. The main advantage of using the definitions (IAip - (IA3p is that if we write k in 
spherical polar coordinates {k, 9, 0) as 

k = k{ei sin 6 cos 0-1-62 sin 6 sin (p + cos6), (All) 

then it is easy to see that 

Ci (k) = Ci cos 6 cos 0-1-62 COS 6 sin — 63 sin 6, (A12) 

e2(k) = — 6i sin0 -|- 62 COS0, (A13) 

63 (k) = 61 sin 9 cos 0-1-62 sin 6' sin -|- 63 cos 9. (A14) 

Now, from Eqs. (1A12|) - (]G13P it immediately follows that the basis {6i(k), 62(k), 63(k)} 
coincides with the spherical basis {65), 6^, 6^}, namely 

ei(k)= l-QQ^^o, (A15) 

/, X 1 <9k . . s 

k sm 9 0(j) 
dk 

e3(k) = ^ = e,. (A17) 
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Finally, we stress that the choice of the recipe (lAltlA3|) is determined by the necessity of 
having two unit vectors parallel and orthogonal to the plane of incidence common to k and 



APPENDIX B: FOURIER TRANSFORM AND BASES CHANGE 

In this appendix we shortly review a few elementary facts about Fourier transform and 
basis changes. To begin with, let us study how the Fourier transform F{pi,p2,P3) of a 
given smooth function f{xi,X2, x^) changes, when we pass from the initial set of coordinates 
{xi,X2,X3} to a new one {ui,U2,us}, via a linear transformation. By definition of Fourier 
transform, the two functions f{xi,X2,X3) and F{pi,p2,P3) are related by the formulas 

f{x,,X2,Xs) = j F(pi,p2,P3)e'^"^P^+"^^^+"^^^^dpidp2dp3, (Bl) 

F{pi,P2,P3) = j f{x,,X2,xs)e-'^^^P'+^'P'+^'P'Mxidx2dxs, (B2) 

with Xi eM., Pi E M, {i = 1, 2, 3), and we use the convention that all integrals are calculated 
in the interval (— oo, oo), unless otherwise stated. The validity of the equations above can be 
easily checked by substituting Eq. ( ]B2I) into Eq. ( IBll) . interchanging the order of integration, 
and then using the following integral representation for the Dirac delta function: 

5(3) (x - x') = j e'P<^-^'M^x, (B3) 

where x = {xi,X2,X3), p = (^1,^2,^3), = dxida;2dx3, and the dot " • " denotes the 
ordinary scalar product in M^, e.g., x ■ p = XiPi and, from now on, we sum over repeated 
indices unless otherwise stated. Now, let {ui, ^2,^43} be a new set of coordinates related to 
the old ones {xi,X2,X3} via a non-homogeneous orthogonal transformation: 

Xi = DijUj + Qi X = Du + a, (B4) 

where, by hypothesis, DD^ = I. Under this transformation, the infinitesimal volume x 
becomes 

dxi dx2 dx3 = (det J) dui du2 du^, (B5) 
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where J : Jij = dxi/duj = Dij, = 1,2,3), denotes the Jacobian of the transformation. 
If we use Eq. (IB4|) in both sides of Eq. (IB II) . we easily obtain 



(2vr)3/2 J ^^P^^ 



(27r)3/2 

where | det J| = | det D\ = 1, and 



j [F{D]i)e'^-^^]e'''-^d^k, (B6) 



■ p = (-Du)iPi = DijUjPi = UjDjiPi = u ■ D^p. 



Moreover, in the third hne of Eq. (]B6P we have changed integration variables passing from 
p to k : k = D^p with, once again, | det J\ = \ det D\ = 1. On the other hand, if we define 

f{Du + Si)=g{u), (B7) 

then, by definition of Fourier transform as given in Eq. fIBip . we can write 



^ j G(k)e^-M3fc. (B8) 



' (27r) 

At this point, it is enough to equate the right sides of Eqs. (IBGHBSp to find the sought 
relation: 

G(k) = F{Dk)e'''-^^. (B9) 

APPENDIX C: GAUSSIAN BEAMS AND BASES CHANGE 

Let K = {Oxyz) be a given Cartesian reference frame, and let Tq denotes a given vector 
pointing along the direction Vq : Tq = rgro, where Tq = |ro| and |ro| = 1. Moreover, let 
{ci, 62, 63} = {e^, By, Bz} be three mutually orthogonal unit vectors pointing along the axes 
x,y and z, respectively. The two vectors Vq and 63 determines a plane, say, the plane of 
incidence. If we think of kq = —r^ as a unit vector pointing along the axis z' of a new frame 
K' = {O'x'y'z') centered in fq, as shown in Fig. CI, we are free to choose the directions of 
the axes x' and y' around the axis z'. By convention, we choose x' on the plane of incidence 
and y' orthogonal to it. Then, if we denote with {e[,e'2,e'^} three unit vectors pointing 
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along the axes x', y' and z\ respectively, the convention above is automatically satisfied by 

choosing, according to the results of Appendix A, {e'l, 63, 63} = {ei(/«o), e2(/«o), 

Now, let P an arbitrary point of coordinates OP = r = (xi, 0:2, 3:3) = r in X, and 
> 

coordinates O' P = r' = (x'^, Xg, Xg) = x' in K' , where the symbol "=" stands for: "is 

represented by". Moreover, the coordinates of the origin O' of K' ^ with respect to K, 
> 

are equal to 00' = ro = (a^oi, 2:02, 3:03) = i"o, by definition. Note that, in order to avoid 
confusion, we use the symbols r, Tq to denote both the vectors and their representation in 
the frame while we use the symbol x' for the representation of the vector r' in the frame 
K'. 
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(/ 
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' ■iiisiii 

iW||||B 





X 



FIG. 3: Coordinate system definition. The two vectors e'2 and ei^ are directed towards the reader, 
parallel to By. 



Since, trivially, 

00' + 0^ ^ OP ^ ro + r' = r, (CI) 
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then, we can write 



= e' ■ r' 



e[-ir- ro) (C2) 



3 



3 



where in the last hne we have defined Dij = e[ ■ ej. If we define 

tti = -DijXoj = -e'- ■ ro = roe[ ■ kq = ro^^s, (C4) 
then we can rewrite in compact form 

X- = DijXj + ai = e[-r + ro^js, (C5) 

or, equivalent ly, 

x' = D{r-vo) = Dv + a, (C6) 

where, in analogy with Eq. (lB4p . we have defined a = —Dtq. It is easy to prove that D is 
orthogonal, that is DD^ = /: 

{DD^)ij = D,kD,k = [e'i ■ efc)(e;. ■ e^) 

= 5., = (C7) 

where in the second line we have used the definition of coordinates of the vectors with 
respect to the basis 6^: (eQ^ = ek ■ e[ = e[ ■ e^. 

Now, let us consider a Gaussian beam propagating along the positive direction of the axis 
z' , at wavelength Aq with a central waist wq = w{z' = 0) located at tq. Its functional shape 
is given by 

ij{x',y',z',t)=exp[tko{z' -ct)]f{x',y',z'), (C8) 
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where ko = 27r/Ao and 



f{x',y',z') 



z' -iL 

is a solution of the paraxial wave equation 



exp 



/2 I „/2\ 1 



2 V z'-iL 



(C9) 



+ — + 2iko— ) f{x,y,z) = 0, 



(CIO) 



^ dx"^ dy"^ dz ^ 

with L = koWQ/2 (the Raleigh range of the beam). Let us calculate the Fourier transform 
of fix', y\ z') as 

1 



(27r)3/2 
(27r)V ^L 



fix ,y ,z)e 



P'Mx'dy'dz' 



exp 



2ko 



{pI+pD 



S{P3 + 



pI + pI 

2ko 



(Cll) 



where we have used the elementary Gaussian formula 

/oo 
exp [(ai + ibi) x + (02 + ^^2) a^^] dx 
-00 



7r 



-02 - ^62 



1/2 



exp 



(ai + i6i) 
4(a2 + ib2 



(a2 < 0), (C12) 



with {cij,&i}f=i G M, to calculate 

J(m, y) = / exp ( tux + t 



.ko 
2" 2' - 



1/2 



exp 



2kr 



dx 



[z' - iL) 



(C13) 



Next, what we want to do is to express /(x', y', z') = /(x') in the frame K, that is we want 
to find the function g[r) such that g{r) = / (x' = Dr — Dtq). From Eq. (lC2p it immediately 
follows that 



^(r) 



• (r - ro) - iL 
^ [,fco [e;-(r-ro)]^+[e- -(r-ro)]^ 
I 2 e'3-(r-ro)-zL 

-iL 



e'o ■ r + ro — iL 



X exp 



.fcp (e;-r)^ + (e^2-r 
2 63 ■ r + ro — iL 
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(CM) 
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and that 



exp[iko{z' — ct)] = exp[ikoe^ ■ (r — fq) — ikoct] 
= exp[iko • (r -ro) - iujot], 



(C15) 



where, by definition, ko = koe'^, and ujq = k^c. Thus, we can rewrite our initial Gaussian 
beam ilj{'x',t) in the frame K as ilj{x.',t) = ^(r,t), where 



^^(r, t) = g{r) exp[iko ■ (r - fq) - Uot]. 



(C16) 



It is well know that this field is not a solution of the Maxwell wave equation, that is 

Q2 Q2 Q2 I d'^\ 

+ — + — - ^(r,t) = □^(r,t) ^ 0. 



dx^ dy"^ dz^ dt"^ ^ 
However, we can seek a field \E'o(i",t) such that □\l/o(i",^) = 0, Vt > 0, and 

*o(r,t = 0) = ^(r,t = 0). 



(C17) 



(C18) 



In order to determine \I^o(r, t), we begin by writing it as a plane wave expansion of the form 

1 



^o(r,t) 



(27r)3/2 



a(k) exp[i(k • r — a;t)]d /c 



(C19) 



where uo = |k|c, and the amplitudes a(k) have to be found. Then, we perform the following 
change of variables: k — * q, where 



k = ko + (k - ko) = ko + q, d^fc = d^g. 



(C20) 



Thus, Eq. fIClQP becomes 



^o(r, t) = j «(ko + q) exp[i(ko + q) ■ r - iuj{\iQ + q)t]d^g 



(27r)3/2 



a(ko + q) exp[iq • r — i5co'(ko, q)t]d'^g. 



(C21) 



where we have defined 



5(x;(ko, q) = c^(ko + q) - c^o 



2\ 1/2 

cfcol 1 + 2^.0-^ + 1 



a;o < kq ■ 



q 



k"^ 
^0 



Ko ■ 



q 



1 - Ko ■ 



q 



(C22) 
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with = ko/fco = 63, and q = |q|. The approximate equahty in the last hne of Eq. (1C22P 
holds in the limit |q|/fco -C 1- At t = Eq. flC2ip can be written as 

„ikor 



*o(r,t = 0) 



a(ko + q) exp(iq ■ r)d g. 



(C23) 



(27r)3/2 

By remembering Eqs. (IClSp and (lC16p . and simplifying the common term exp(iko ■ r) on 
both sides, we can rewrite Eq. ( ]C23p as 

1 



g{r)exp{-iko ■ fq) 



(27r)3/2 

which, according to Eqs. (lBltlB2p . implies 



a(ko + q) exp(zq ■ r)d g. 



a(ko + q) = G'(q) exp(-iko ■ Fq). 



(C24) 



(C25) 



were G{q) denotes the Fourier transform of g{r). Now, to complete the calculation we have 
to remember that we have defined g{r) = f{Dr — Dtq) and that, according to Eq. (lASp we 
can write 



G(q) = F(Dq)exp(-^q-ro), 



(C26) 



where Eqs. flB4|C4p have been used. By substituting Eq. flC26P into Eq. flC25P we obtain 



a(ko + q) = F{Dq) exp [-i{ko + q) ■ fq 



(C27) 



This equation can be now used in Eq. (]C2ip to obtain our final result 



^o(r,t) 



gi(ko-r-ajot) 

(27r)3/2 
(27r)3/2 

gi[ko-(r-ro)-u;ot] 

(2^0^72 
gi[ko-(r-ro)-iJO<] 



(27r)3/2 

where in the last line we made the change of variables 

q^p = Dq. 



a(ko + q) exp[iq ■ r — i(5u;(ko, q)t]d'^g 
F(Dq)e-^(''o+'>)-'"o exp[iq • r - i6uj{ko, q)t]d\ 
F{Dq) exp[iq • (r — Fq) — i6Lj(ko, q)t]d'^q 



From Eqs. (EI]), (IC2]), and (IC22!l . it readily follows that 

{D^p) ■ (r - ro) = p • D{r - fq) = p ■ x', 



(C28) 



(C29) 



(C30) 
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and 



UJQ 



1 + 2ko ■ 



1 + 2^ + ^ 
ko kl 



ko V ^0 

2\ 1/2 



1/2 



UJo 



1 



where {D'^pf = (D'^p) ■ (D'^p) = pDD^p = pp=p^, and 

Ko ■ (D'^p) = (Dko) ■ p = Dij{e'^)jPi = DijDj^pi = p^. 



(C31) 



(C32) 



Since the Dirac delta function in Eq. (IClip imphes that 2p^/ko = (p| —p'^)/kQ, it follows 
that we can rewrite 600 {D^p) as 

.2 \ 1/2 



SLj{ko,D p) = uo 



LOo 



UJo 



ko k^ 



„2x 1/2 



/ 2 I 2\ 2 
^0 f P 1+P2 

2^0 



+ 



(C33) 



where the last approximate equality holds in the paraxial limit |p|/fco 1- Thus, we see 
that for an input Gaussian beam, the non-paraxial corrections to the carrying frequency uJo, 
begin at fourth order. Finally, by using Eq. (1C30I) . we can rewrite Eq. (1C28I) as 



gi[ko-(r-ro)-wot] 

*o(r, t) = jT^z^, / F{q) exp [iq ■ x' - ir(q)t]d=^g. 



(C34) 



(27r)3/2 

where we have defined r(q) = 6uj(ko, D'^q). It is gratifying to see that after this tour de 
force we re-obtained our starting result (1018^ : 

iko-{r-ro) 



*o(r,t = 0) 



F(q) exp {iq ■ x )d g 



(27r)3/2 

giko-(r-ro) j^^/ = - Dro) 



^(r,t = 0). 



(C35) 
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1. Geometric reflected scalar beams 



Now that we learned how to represent a scalar Gaussian beam in an arbitrary Cartesian 
frame, we can study how such representation changes under reflection. In Sec. I, we showed 
that a generic scalar fleld 0(r) transforms to 0(r) = 0(r) under geometric reflection with 
respect to the plane of equation ^ = 0. Therefore, by using Eq. (]C34I) we can write at once 
^o(r,t) ^o(r,t), where 

_ gi[ko-(r-ro)-wot] f 

^o(r, t) = j^-^^ J F(q) exp [^q ■ (Dr - Dvo) - tT{(i)t]d% (C36) 

However, in this form, this equation is of little usefulness. In order to gain insights, flrst we 

a) deflne a "mirror image" local frame K" specular to K' (as shown in Fig. 3), and then 

b) we see how a scalar fleld represented by (j){x,y,z) in K, and by (p' {x' , y' , z') in K', can 
be represented by (j)"{x" ,y" , z") in K", where (f)" {x" , y" , z") is to be found. The Cartesian 
frames deflnitions are illustrated in Figs. CI and C2. 

Thus, let K" = {0"x"y" z") be a Cartesian frame centered in ro = — ro/5o, and let 

{e'/, e'g, e'g} = {ei(K.o), e2(/5o), e3(K.o)} be three unit vectors pointing along the positive 

directions of the axes x", y" , and z" , respectively. They are deflned via Eqs. (]AltlA3P 

with the substitution k kq. As before, let P denotes an arbitrary point of coordinates 

OP = r = (xi,X2,X3) = r in i^, and coordinates 0"P = r" = {xi,X2,x'^) = x" in K". 

> ^ 

The coordinates of the origin O" of K", with respect to K, are equal to 00" = fq = 
(xoi,Xo2,Xo3) = (a;oi,a;o2, -xos) = r o where, by deflnition, fq = (xqi, Xo2, Xqs). 



Then, the following equality it is trivially true 



00" + 0"P = OP ^ vo + r" = v. (C37) 
Therefore, we can always write 



x" = e" ■ r" 



= <-(r-?o) (C38) 

3 
3 

= ^Dij{xj -xoj), (C39) 
i=i 
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FIG. 4: Illustrating the definitions of some vectors used in the text. 



or, in vector notation 



x" = D(r-?o). 



(C40) 



where we have defined Dij = e" ■ ej. It is easy to prove that D is orthogonal, that is 
DD^ — I. Now, without loss of generality we choose the Cartesian frame [Oxyz) in such a 
way that both Kq and lie in the xz plane (the plane of incidence), that is we write 



Kq — Ci sin ^0 + 63 cos 
Ko = ei sin ^0 — ^3 cos ^o- 



(C41) 
(C42) 



In such a frame it is straightforward to calculate: 



D 



( cos — sin ^ 
1 
sin^o cos^o 



(C43) 



D 



^ — cos ^0 — sin ^0 ^ 
1 
sin^o — cos^o 



XDZ, 



(C44) 
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where we have defined the three "refiection" matrices X, Y, Z as 



X = 



\ 



-10 
1 
1 



, Y = 



1 
0-10 
1 



, z = 



1 
1 
0-1 



(C45) 



Since = = = I, from Eq. fICiil) it readily follows that D = XDZ. Then, given 
any vector u = {ui, U2, Ms) in K, one can use Z to write 



u = Zu, 



from which it trivially follows that 



Du = DZu = X'^DZvL = XDu, 
Du = X^DZ'^u = X^DZu = XDu. 



(C46) 



(C47) 
(C48) 



Note that although we have derived the relations flC44IIC47|C48p starting from the definitions 
(lC4mC42l) . they are independent from the latter. At this point we can go back to Eq. (1C36P 
and notice that we can rewrite 



and 



ko ■ (r - ro) = ko • Z{r - fq) = Zko ■ (r - fq) = ko ■ (r - fq), 



Dr - Dro = XDr - XDvq = XD{y - ?o) = Xx", 



(C49) 



(C50) 



where Eq. 
obtain 



has been used. If we substitute Eq. 



into Eq. flC36p . we finally 



^o(r,t) = 77::z^9 / -P'(q) exp [iq ■ X£)(r - Fq) - ir(q)t]d^g 



(27r)3/2 

gi[ko-(r-ro)-a;oi] 



F(q) exp [iq ■ Xx" - ir(q)t]d^g. 



(27r)3/2 

It should be noticed that at t = such expression reduces to 

$o(r,t = 0) = e^'=o^7(Xx"), 
where we used Eq. ( ]C38I) to write 



(C51) 



(C52) 



ko ■ (r - ro) = huo ■ (r - ro) = hel ■ (r - ro) = kox^ = koz" 



(C53) 
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Equation (1C52|) is exactly what we expect for a Gaussian beam centered at r^o and directed 
along ko- It is also important to notice the argument Xx" of the function / that cor- 
rectly account for the parity inversion between local frames in reflection. In other words, if: 
(f){x, y, z) and y, z) = y, —z) are a field and its mirror image as seen from the frame 
K\ y', z') is the field z) as seen from the local frame K' , and, finally, (f)"{x", y", z") 

is the field (f){x,y,z) as seen from the local frame K"; then 4>" {C, , tj , () = 4''{~^,vX)> where 
(^, 7], () is a set of three given real numbers: 

0(r) = 4>{Zv) 



K' 



K" 

'(Xx') 



where x' = Z)(r — fq) and x" = D{r — tq) 



APPENDIX D: VECTOR BEAMS 



(C54) 



Up to now, we considered only scalar beams of the form 

1 



^(r,t) 



a(k) exp[2(k ■ r — ujt)]d k, 



(Dl) 



(27r)3/2 

which could be a faithful representation of an unpolarized beam. However, if one is interested 
to polarization properties of the field, such representation does not work and one has to pass 
to a vector description. Since a polarized beam can be physically obtained on a laboratory 
bench by sending an unpolarized beam across a polarizer oriented, for example, at /, one 
could naively think that it would be possible to promote the scalar field \&(r,t) to the rank 
of a vector field ^{r,t), just by multiplying \E'(r,t) by /: 



^(r,t)^*(r,t) = /*(r,t). 



:d2) 



where / is a real- or complex-valued unit vector. This equation is obviously wrong since the 
so-obtained field ^{r,t) is not transverse: 



div*(r,t) 



V • *(r,t) 

' d 



4 = 1 



(27r)3/2 
7^ (in general). 



j a(k)(k ■ /) exp[z(k ■ r - Ljt)]d^k 



(D3) 
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However, for any wave vector k we can always write / as 

/kk\ , / kk\ , 

= /:(k)/ + T(k)/ 

^ f||(k) + f^(k), (D4) 

where f±(k) is genuinely transverse, namely k-f_|_(k) = 0, and Eq. flASp has been used. This 
fact suggests us the possibihty to substitute each scalar amplitude a(k) in the plane-wave 
expansion (IDip . with the vector amplitude a(k)fx(k): 

a(k) a(k)fx(k), (D5) 

thus obtaining a bona fide transverse field 

*(r,t) = J a(k)f^(k)exp[^(k-r-cut)]d3fc. (D6) 

Let {ei(k), e2(k), e3(k)} be an orthogonal and complete basis built following the recipe 
given in Appendix A. By using such a vector basis, it is easy to see that we can rewrite Eq. 
(|D6|) as 

2 

^(r,t) = j^J2j exik)a,ik)exp[t{k.r-ujt)]d% (D7) 
where we have defined 

aA(k)=a(k)e,(k)./. (D8) 

In conclusion, if we have a beam passing through a polarizer oriented at /, it can be described 
by a vector field ^{r,t) such that 

^{v,t) = M{r,t)f, (D9) 

where A^(r, t) is a 3 x 3 matrix defined as 

Mir, t) = j T(k)a(k) exp[z(k ■ r - ujt)]d^k. (DIO) 

Now, we apply the previous formulas to the actual case of a quasi- mono chromatic Gaus- 
sian beam propagating along the positive direction of the 2;-axis, whose equation is 

V^(r,t) =exp[^A;o(^-ct)]/(r), (Dll) 
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where 



/(r) 



-iL 
z — iL 



exp 







I — 
2 


[z-iL J_ 



(D12) 



Here, Aq = 2TT/ko is the carrying wavelength of the beam, c is the speed of hght in vacuum, 
and L = koWQ/2 is the Raleigh range of the beam whose waist at z = is equal to: 
w{z = 0) = wq. In this Appendix we are not interested to the temporal evolution of the 
beam, so we assume t = in all subsequent formulas. Then, we rewrite Eq. (IDip as 



^(r,0) 



(27r)3/2 



a(k) exp[ir ■ (k — koez)]d^k 
a(q + koCz) exp(2r ■ q)d^g. 



(D13) 



(27r)3/2 

where in the second line we have changed the variables of integration from k to q = k — fcoe^. 
If we impose the condition \l/(r, 0) = tp{r, 0), then from Eq. fIDllI) and Eq. fIDlSI) . it readily 
follows that 



/(r) 



a(q + koBz) exp(zr • q)d^q, 



(D14) 



(27r)3/2 

which implies that a(q + kQe^) coincides with the Fourier transform F{q) of /(r): a(q + 
/coe^) = F(q), namely 



a(q + koB^ 



(27r)V2L 



kn 



exp 



2ko 



5 U3 + 



ql + ql 

2kn 



(D15) 



where Eq. fIClip has been used. Note that the presence of the Dirac delta in the Fourier 
transform -F(q) of /(r) is not accidental. In fact, given a function /(r) solution of the 
paraxial wave equation: 



— + — + 2tk — 



f{r) = 0, 



(D16) 



it is always possible to write 
= 



q2 q2 

9x2 + + ^ "^dzj [(27r)3/2 

-2kn 



F(q) exp(ir • q)d^g 



(27r)3/2 



F{q) qs + 



qj + gj 

2k. 



exp(?r ■ q)d^g. 



(D17) 
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This equation implies that 



2 2 ^ 

i^(q)(g3 + ^i^) = V(gi,g2,g3)e 



(D18) 



namely F{q) must be zero outside the 2D surface of equation + ^^^^ = 0. However, in 
order to have /(r) ^ 0, it is clear that F{q) must be infinite upon such surface. Thus, it 
must be always possible to write 



F(q) = J-(q) S(qs + 



oi + ql 



2h 







(D19) 



where J^(q) is finite everywhere in R^. At this point it is a trivial task to rewrite Eq. (IDlOp 
as 



-^('''O) = j ^(q + ^oej;)a(q + fcoe^) exp(ir ■ q)d^g 

dgidg2S exp 



(27r)3/2 



27r/Cn 



ia;)gi + {iy)q2 - T^iqf + ql) 

ZKq 



X J T(gi, ga, gs + ^o) exp (izq^) 6 [ qs + 



2TTkn 



exp 



{ix)qi + {iy)q2 



ql + ql 

2kn 



L -\- iz . 2 2\ 



(gi, 92, /i;o - (gi + qi)/{2ko)) dgidga, 



(D20) 



where the symmetric transverse tensor T(q + fcoe^), after integration with respect to the 
variable q^, takes the value 



T{quq2,ko-{ql + ql)/i2ko)) 



1 



_QiQ2_ 

l+Ql/4 



Q1Q2 



l+Ql/4 



2i(l-Qi/2) \ 

l+Ql/4 
?2(l-Qi/2) 



(l-Qi/2) _Qi(l-Qi/2) 
\ l+Ql/4 



l+Ql/4 
n2 



( 



-Q1Q2 -Q 



I+QYI^ l+Ql/4 / 



(D21) 



-Q1Q2 -g2 

\ -Ql -Q2 Ql + Ql) 

In the first line of the equation above, we have defined the dimensionless variables {Q\ = 
Q'a/^o}a=1) Q\ = Ql + Q\\ while in the second line we have made a Taylor expansion 
about Ql = = (^2 up to the second order, which is justified for paraxial beams. 
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To illustrate in detail the formalism just developed, let us apply it to calculate the vector 
structure of a Gaussian beam passing trough a polarizer oriented along the x-axis: f = e^. 
From Eq. it follows that 

*(r,0) = M{r,0)e, 



1=1 



Mu{r, 0)e^ + 7W2i(r, 0)ey + Msiir, 0)e, 



(D22) 



The matrix elements Aiij can be easily calculated in the second order approximation by 
using Eqs. flD20tlD2ip as follows: 

A^n(r,0) ^ / exp[itkox)Qi + itkoy)Q2] {l - Ql) 



271 



1 



xexp <j [koL + (tkoz)] {Ql + Ql) )■ dQidQa 



271 

xexp 

271 

xexp 



J exp (XiQi+Xags) {l - Ql) 

dQidQ2 

(1 - Ql) 



--(A + X3) {Ql + QD 



1 ^ _ 2 



g,/3=l 



a=l 



d^Q, 



(D23) 



where we have defined the dimensionless variables 



Xi = ikox, X2 = ikoy, X3 = ikoz, A = koL, 



and the covariance matrix A of elements 

1 



A. 



{a,Pe{l,2}). 



(D24) 



(D25) 



(A + X3) 

Gaussian integrals are easy to calculate and the following result can be found in any graduate 
textbook: 

1 ^ _ ^ 

■■^ ^ QaiA ^)a/3Q/3 + QaXg 



G 



exp 



a=l 



d'Q 



27r(detv4)^/2 exp 



I 2 

9 ^g^g/3-^/3 



g,/3=l 



(D26) 
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In our case, with A given by Eq. (lD25p . we obtain 



G 



27r 



exp 



1 f xf + xl 

2 V A + Xs 



(D27) 



By using this result we can easily calculate the following two quantities that we need: 



Gn 



f r 1 ^ - ^ 

L a,/3=l a=l 

dG 



G 



and 



Gad 



A + Xs' 

(5q(5/3 exp 
d'^G 



(D28) 



1 ^ _ 2 



9X„9X^ 

x^x^ 



+ 



a,/3=l 



a=l 



d^Q 



(D29) 



[(A + X3)2 A + X3J 

Now we have all the ingredients we need to calculate the three matrix elements required to 
evaluate Eq. ( ]D22I) . In sequence, first we calculate A^ii(r,0): 



A^ii(r,0) 



2n 
-iL 
z — iL 



G 



exp 



+ . . . 



d^G 
dXf 

ikoz + 1— 



X 



1 + 



X 



z — iL 
2 



^(r,0) 



z — iL \z — iL 

1+' ' 



+ 



kn z — iL \ z — iL 



+ 



(D30) 
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where Eq. ( ]DimD12|) have been used. Then we calculate A^2i(r,0): 



A^2i(r,0) 



2ti 
xexp 



exp (XiQi + X2Q2) (-Q1Q2) 
4(A + X3) {Q\ + Ql)\dQ,dQ2 



z — iL 
^(r,0) 



2n 8X28X1 
iL 

exp 



•7 , -^0 

ikoz + 1— 



z — iL 



xy 



[z - iLf 



+ 



xy 



[z~iLf 



+ ... 



and, finally, we calculate A^3i(r,0): 



A^3i(r,0) = ^ / exp(Xigi+X2Q2)(-Qi; 



27r 
xexp 



(A + X3) (Q? + g^) UQidQ2 



271 
iL 

exp 



z — iL 
V'(r,0) 



■ .fco f x^ + y^ 

ikoz + I— — 

2 \ z — iL 



X 



z — iL 



+ 



a; 



(D31) 



(D32) 



z — iL 

By collecting all the results above we can obtain the final expression for a unpolarized 
Gaussian beam crossing a polarizer oriented at f = Bx, exact up to the second order terms: 

2" 



*(r,0) ^^(r,0) <^ e. 



1 



1 



ko z — iL \z — iL 



xy 



^'{z-iLf 



X 



z — iL 



. (D33) 



As expected, the dominant zero-order term is the "naive-guess" vector field ^/^(r, 0)63,. The 
first order correction gives a longitudinal contribution oc to the field. Finally, the so- 
called "crossed-polarization" term oc amounts to a second order correction. A plot of the 
intensities of the three terms is reported in Figs. 5-7. These results should be compared 
with the similar ones presented in ref.— . 
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FIG. 5: Transverse intensity profile due to the first term oc Sx (dominant polarization term) in Eq. 




FIG. 6: Transverse intensity profile due to the second term oc (cross polarization term) in Eq. 
(ID33]) . 




FIG. 7: Transverse intensity profile due to the third term oc (longitudinal polarization term) in 
Eq. (ID331) . 
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APPENDIX E: EXPANSION OF THE POLARIZATION VECTORS 



Let {ei(k), e2(k), e3(k)} be an orthogonal and complete basis for defined as in Ap- 
pendix A, where k G is a given vector. Each basis element ej(k) can be expanded in a 
Taylor series around an arbitrary vector k = kg. If kg denotes the "carrying" wave vector 
of a well collimated light beam, then it is reasonable to keep only the first and second order 
terms in the Taylor expansion, and write 

5e,(k) 



ei(ko + q) ~ ei(ko) + qi 



dh 



k=ko 



^ 1 g^e,(k) 



+ 



k=ko 



,(0) 



.(0 



[Im) 



+ 



where we have defined 



?^-ef)(ko 



ef^er^(ko) = fco 



ei(ko), 
9e,(k) 



dh 



(Im) 



(Im) 



(ko 



k=ko 

2 d'e^ik) 



(El) 

(E2) 
(E3) 

(E4) 



k=k 



■0 



'° dkidk 

Since we have always the freedom to choose the right-handed basis {e^, Sy, Sz} = {ei, 62, 63} 
that defines orientation of the Cartesian frame (Oxyz), we can take ko laying on the xz- 
plane: 

ko = ko (e^ sin 9o + cos 9q) = ko 



(E5) 





^cos 6*0^ 

where ko = |ko|. In this basis is a straightforward task to calculate, by using the definitions 
given in Eqs. (]E2tlE4p . the three zero-order vectors {e[^^}f^i, the nine first order vectors 
{ef^}f and the twenty-seven second-order vectors {ef™^}f ; j^^^. Explicitly, the zero-order 
vectors are 



/ .... \ 



COS 6^0 


y— sin 6*0 
The first-order vectors are 

^ cos 6*0 sin 6*0^ 



^2 





COS^ 00 



1 







sin 9o 


ycos 6*0 y 



(E6) 



COS^ 60 



cos 6*0 sin Oq 



(E7) 
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/ \ 



ef = 





cot 

V I 

( sin^^o ^ 


cos ( 



,(2) _ 





V ° / 






1 





(E8) 



^— cos ^0 sin ^o"^ 



Since e|'"^^ = ej'"'^ , the o: 

^ ~\ (cos ^0 + 3 cos 3^o) 


^ 3 cos^ 9q sin j 







(E9) 



sin^ 6*0 J 



ycos ^0 sin 6*0 J 

)nly relevant second-order vectors are 

/n\ / -3cos2^osin^o 



„(11) _ 



p(12) _ 



^ N 

— COS 9o — cot 00 esc 00 




V 



e(12) _ 



p(13) _ 



(1 + 3cos2^o) sin^o 


I (cos ^0 + 3 sin 36*0) j 



p(13) _ 







w 

/csc2^o^ 



V J 

,(13) 







^— I (cos ^0 + 3 cos S0o)j 



p(22) _ 



\4 
/ 

— COS 6*0 — cot 6*0 CSC 0o 


^ — cos 00 cot 00 J 

( \ 



5 ^3 



V7 



63^^'' = - sin 00 

\ " / 

^— \ (cos ^0 + 3 cos 3^0) 


— I (sin 6*0 — 3 sin 36*0) y 



(ElO) 

, (Ell) 

A 



V 



p(22) 





— CSC 





■'00 



/-sin 



(E12) 



*^3 — 



V ^ / 



p(23) _ 



p(33) 



^ —3 cos 6*0 sin^ 0o 



sin 0Q 

\ 



p(23) 
*^2 







(sin6lo - 3sin36lo)y 



p(33) 
) *=2 








p(23) 
*^3 



— COS ^0 

V ; 



, (E13) 
j^-cos 00 J 

\ 



p(33) 
*^3 



^— I (sin^o ~ 3 sin 35*0) 


— 3cos^o sin^ ^0 j 



(E14) 



A 



(E15) 

Prom the equations above it is very simple to obtain the corresponding formulas for the 
of the polarization vectors defined around the "reflected" vector ko, by noting 



expansions 
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that 



and 



and 



ef . 






62 = 


.e<»)(k„) 


- 

~ "=2 ' 

- - 6^°) 
~ ^3 




- e»(k„; 




^3 


= ^2 IKo^ 


— ^2 , 


(hn) 


e^)(ko) 


= - ei 


(hn) 

2 — 

(hn) 

3 — 


er(ko) 


~ "=2 ' 

= - 4 



6*0 -^-9o 

{]m\ (]m\ \ (]m\ 

e 

A quantity that often appears during calculations, is the diadic product 



(E16) 

(E17) 
(E18) 

(E19) 

(E20) 
(E21) 

(E22) 

(E23) 
(E24) 



e^(ko + q)e,(ko + q), (/i, z/ G {1, 2}), (E25) 
which can be expressed in terms of the expansions ( lEip . up to second order terms, as 



e^(ko + q)e,(ko + q)~eWeW + 



(O)p(O) , ^ 



p{0)p(i) 



+ 



2^0 



mi 



e-e--+2(-ir«eWei^^ + (-l) 



.(O)pfe') 
1 



-J-a+^jSpfoOpCO) 



E^, + g,E« + -Q,Q,E(:;^), 



(E26) 



were we have introduced the dimensionless variables {Qi = gi/fco}f=i- The 4 + 12 + 36 
matrices (each 3x3) E^^, E^^*,^ and e|i*^\ can be calculated straightforwardly from the 
following definitions: 





= p(0)p(0) 








= p(0)p(ii) 



'5«o(i)p{0) 



-l)""eh'''e 



r_l)'5»3+<5i3efe')e(o)_ 



(E27) 
(E28) 
(E29) 
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Note that the presence of terms hke cot Oq and esc 60 make the expansion above not 
utihzable for "small" values of 60 . More precisely, it is possible to show that the expansion 
(lEip breaks down when 

• n ^ ■ /)_ |ko X (ko + q)| , - 

sm6'o<smt^ = — — 1 — , (E30) 

^o|ko + q| 

where Eq. (2.9) of ref.— has been used. 



APPENDIX F: PARAXIAL FORMULATION 



In this appendix we consider an alternative and simpler formulation for the study of 
the reflection of a Maxwell-Gaussian beam by a planar surface. Basically, we use here the 
formalism developed by Haus and Pan in Ref.^^. The geometry of the problem and the 
notation we use are those ones fixed in appendix C. Let A™(x', y', z', t) the vector potential 
of the impinging electromagnetic field, as seen from the local frame K' . By hypothesis, it 
obeys the Lorentz gauge condition. We assume that A™(x', ?/', z', t) can be written in the 
paraxial approximation, as 

A-(a:', y\ t) = [e^ij^i^') + e',^'^{^')] exp{-tu;ot) (Fl) 

where, from Eq. (C6), x' = D(r — fq), and 

C(x') = uJ{x',y',z')exp{zkoz'), (a = 1,2), (F2) 

where f{x', y', z') is defined by Eq. (C9), and G C : + |u2p = 1. Following appendix 
C, we define gir) = fix.' = Dr — Dyq), and use Eqs. (C15,C16,C18) to rewrite 

^lJ':^{x' = Dr - Dro) = u^^iv,0) 
= M„^'o(r,0) 

= Ua (-2^)3/2 J ^('^o + q) exp(iq ■ r)d^g, (F3) 

where Eq. (C28) has been used. 

Let A°^^ [x" , y" , z" , t) the vector potential of the reflected electromagnetic field, as seen 
from the local frame K": 

A°-\x", y", z", t) = [e;Vr (x") + e'ir2^\x")] exp{-^u;ot). (F4) 
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Since reflection causes r — > r and a(ko + q) — > ra(ko + q)a(ko + q), we can write at once 



r^iko + q)a(ko + q) exp(iq • r)d^g, (F5) 



^(27r)3/2 

where there is no summation over repeated Greek indices. At this point, it is easy to follow 
the procedure outlined in Sec. Ill in order to obtain 



U 



(27r)3/2 

ikoz 



j r„(ko + /^^q)F(q) exp(iq ■ Xx")d^g 



s^'^"^ a„(x") 



Zn 



(F6) 



where Ga(x") is deflned by Eq. (53). From Eqs. (8-9) of Ref.— and the equation above, it 
is easy to calculate the time-independent parts of the reflected electric and magnetic flelds: 

2 



Efx") 



ikoz'' 



a=l 
2 



B(x") = zkoe''''^" V (e^; x e:)a^ + - (e'^ x V>«) 



where 



d d 



dx" 



dy" 



a„ X 



(F7) 
(F8) 

(F9) 



These expressions are both analytical and easy to handle. They can be used for the flnal 
step, namely the computation of the time average Poynting vector S(x"): 

1 



S(x") = —Re [E(x") X B*(x" 
2/io 



(FIO) 



The flux of the Poynting vector through the detector's surface gives the measured intensity. 
Assuming that the normal to the detector's surface is directed parallel to 63, we have to 
calculate the position of the center of the function 5*3 (x") at the position z" = of the 
detector: 



\a,ix",y",Zn)\' + \a2ix",y",Zn)\'). 



(Fll) 



Thus, the GH and the IF shifts, as measured in K", are simply: 

2 



j: /.i.„(.,.,z„)iM.d. 

a=l 

2 ' 

^ / \aa{x,y,ZD)\^dxdy 

a=l 



(F12) 
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^ j y\aa{x,y,ZD)\^dxdy 

LiF = ^ , (F13) 

^ / \aa{x,y,ZD)\'^dxdy 

respectively. Note that these integrals are all Gaussian, so that they can be calculated 
analytically. 

It is quite clear that this approximate formulation can give a correct expression for the 
GH shift, but not for the IF shift. 



APPENDIX G: POLARIZATION OF LIGHT BEAMS 



Notation 



The scalar product between two complex-valued three-dimensional vectors u and v is 
defined as 

3 

^' ^ X] '^i^i ^ ^i'^i ^ ('^' u)* eC, (Gl) 



3 

[U, 

i=l 

where summation over repeated indices is understood, and Ui, Vi are the Cartesian component 



of u and v, respectively, with respect to an arbitrary orthogonal reference frame. Unit vectors 
are denoted by the " hat " symbol: 

u^^^^-, {u,u) = l, (G2) 

where u = {u, u) . The diadic product uv'^ represents a 3 x 3 matrix whose elements are 
defined as 

{uv'^)ij = UiV*. (G3) 

Operators are denoted by capital letters either as P or P; their matrix representations in a 
given basis are always written in "bold" characters as P. 



Arbitrary elliptical polarizers 

Let us consider the action of a polarizer upon a quasi-monochromatic well coUimated 
beam of light that crosses it non-orthogonally. Let {Oxyz) be the laboratory frame specified 
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by the real-valued basis 

{x,y,z} = {61,62,63}. (G4) 

An elliptical polarizer on the laboratory bench is characterized by its axis z and its orien- 
tation p: 

P = PxX + Pyy, (G5) 

where Px,Py are complex- valued numbers such that jpxP + = 1- This vector represents 
an optical device made by three optical elements: a quarter wave plate, a linear polarizer, 
and another quarter wave plate (see Ref.— ). This is easily seen by writing explicitly the 
Jones matrices C(a) and P(/?) of a compensator and a linear polarizer, respectively, in the 
two-dimensional laboratory basis {x,y}: 

/ e*" \ [ cos'^p sin/5cos/3\ 

C(«)= . , P(/5)= 2 ■ 

\ 6-'" I \ sin/3cos/3 sin^ (3 J 

These well known matrix representations hold when the impinging light beam propagates 
along the 2;-axis. If, without loss of generality, we rewrite px and Py as 

= e"*" cos (3, Py = 6'" sin /?, (G7) 

then it is easy to check, via a straightforward calculation, the validity of the following 
relations: 

_t I P'-'P*^ P^Py 
pp' = I 

\PyP*x PyP*y 

= C(-a)P(/3)C(a) 

cos^ P 6"^*°" sin P cos P 
e^*° sin P cos P sin^ P 



(G8) 



Now, let's turn back to the three-dimensional problem, and let A; be a real-valued unit 
vector denoting the main direction of propagation of a quasi-monochromatic well collimated 
beam of light whose electric field vector can be written as 

Eo{f, t) = {Ex'X + Ey,y') e^^K^'^^)-^*] = ^e^TK^.-^--*!, (G9) 

where 

{x',y',z'^k}^{6\,6'^,6',}, (GIO) 
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is a local Cartesian basis "attached" to the beam whose axis z' coincides with the beam 
direction of propagation k. If with Q we denote the angles between z and /c, namely (i, A;) = 
cos^, then we can chose the basis {e'j^, e'g, 63} as 

e'l = ei cos Q cos + 62 cos ^ sin — 63 sin ^, (Gil) 

62 = — ei sin0 + 62 cos(^, (G12) 

63 = Ci sin 6* cos + 62 sin 6* sin + 63 cos 6*. (G13) 

— * — # ^ 

The electric field £^ is transverse^ that is (£", A;)=0. In an equivalent manner, we can rewrite 
this condition as an eigenvalue equation 

T{k)E^E, (G14) 

where T{k) is the operator that projects into the space orthogonal to the propagation 
direction k. It is defined as 

T{k)=I-kk. (G15) 

The action of the polarizer upon the field E can be found by requiring that the polarization 
vector T of the transmitted beam lie in the complex plane defined by the polarizer orientation 
p and the propagation vector k of the beam. This condition is automatically satisfied by 
tacking the projection of p upon the space transverse to k, namely 

T — T{k)p — p — k{k,p) = p — ak (G16) 

where a = {k,p) and (r, r) = 1 — lap. Note that r is a unit vector only when the beam 
impinges orthogonally to the polarizer, that is when a — 0. Thus, when the beam of light 
crosses the polarizer oriented at p, its electric field vector becomes parallel to f according to 

E^T{r,E)^TE, (G17) 

where we have defined the operator projector T as 

T = ff^=pp^- {akp^ + a*pk) + \a\^kk. (G18) 

It is easy to check that it satisfies the following properties 

TA; = 0, trT = (f,f) = 1 - (G19) 
T^^T, T^ = (trT)T. (G20) 
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Input and output beams 



So, we have an input beam directed along the axis k whose electric field has Cartesian 
components [E^i^Eyi) in the local "beam" frame {Ox'y'z'). Moreover, the polarizer orienta- 
tion vector p has components {px,Py) in the laboratory frame (Oxyz). Therefore, if we want 
to write Eq. (1G17|) either in the beam or in the laboratory frame, we need to know how 
the electric field and polarizer orientation vector representations changes by passing from a 
frame to the other. Let us define the Cartesian components of the electric field {Ei}^^^ and 
{-Ej'jf^i in the laboratory and in the beam frames, respectively, by the following relations: 

3 3 

E = Y,e^{eiJ) = J2eiEi, (G21) 



i=l 
3 



4 = 1 

3 



E = J2^ie^^E)^J2^El. 



(G22) 



1=1 



i=l 



The change of basis matrix A that permits to pass from one representation to the other is 
defined via its elements as 



or, in explicit form. 



A 



cos 6 cos < 
cos 6 sin <; 
— sin 9 



— sin (f) sin 6 cos ( 
cos sin 6 sin ( 
cos^ 



(G23) 



(G24) 



where AA"!" = A'^A = I. 

— * 

By means of the electric field vector E we can describe any polarized beam of light. 
However, if we need to describe partially polarized light, we must use the coherency operator 
J whose elements are defined in the laboratory and in the beam frames as 

J., = (e„Je,) ={E^E*), (G25) 

Jl^ = {e[Je'^ = {ElE^*), (G26) 

respectively, where angular brackets (■ ■ ■ ) denotes either temporal or spatial (or both), 
averages. From a straightforward calculation it follows that 



J' = A^JA, 
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(G27) 



where, explicitly, for the input beam 

/ 



{Ey'E;,} {e,,e;,} 



V 















(G28) 



Now, following the same line of reasoning, we write the projection operator T in both the 
laboratory and in the beam frames, as 



% = (e:,Te;) = (e:,ffte;) = {e[,r) {r,e'^) ^ t[t^\ 
respectively, where we have defined the representation vectors t,t' as 



(G29) 
(G30) 



f=^eA = ^e:i;. (G31) 

1=1 i=l 

It is easy to see that we can pass from the beam to the laboratory representation of the 
polarizer via the following relation: 



T' = A^TA. 



(G32) 



Now we have all the ingredients to write the relation between input and output beams 
in the frame we want. In particular, from a straightforward calculation it follows that if we 
write such relation as 

Jout = TJinT^ (G33) 
then, in the laboratory and in the beam frames we have 



Jout — TJjnT^, 

t/ rp/ T/ rp/t 

*^out ~ -•- ''in-'- ) 



(G34) 

(G35) 



respectively 



Example 



As an application of the formalism developed above, let us consider the case of an un- 
polarized input light beam directed along k, and represented in the beam frame by the 
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following diagonal coherency matrix 



V 



1 
1 




(G36) 



This matrix is real-valued, so it does not contain any complex phase term. Moreover, let 
suppose to have a circular polarizer represented in the laboratory frame by the orientation 
vector 



(G37) 



By applying Eq. (IGSp we see immediately that in the laboratory {Oxy) frame 

... 1 I I -i 



pp' 



i 1 



C(-7r/4)P(7r/4)C(7r/4). 



(G38) 



Now, from Eqs. ( ]G16IIG24|G31I) it is easy to calculate both t and t , obtaining 
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These equation are interesting for several reasons: First of all, apart from an irrelevant 
phase factor, at zero-order t and ? coincides with a circularly polarized unit vector. The 
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phase factor present in the expression of ? is simply due to the non-univocal definition of 
the spherical coordinate (f) when ^ = 0. Conventionally, in this case one puts = 0, thus 
removing the factor e**^. Second, while in the beam frame the corrections start from the 
second order term in the Taylor expansion of t' around ^ = 0, in the laboratory frame a first 
order (longitudinal) correction term arise. Third, for 6' 7^ 0, the output beam polarization 
vector r has no longer unit modulus: 

(G41) 

Finally, by using Eqs. (!GMT30llG3llG35]lG^^ we can write, after a straightforward 
calculation, the polarization state of the beam seen from both the laboratory and the beam 
frames, as 
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(f, r) = (f, t) = (f, t') = i [3 + cos(2^^)] = 1 - y + 0{df. 
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2cos(20) f i - 2 sin(2(/)) 



fz-2sin(20) -f + 2cos(20) 
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As expected, both representations give a pure circularly polarized beam at zero order (re- 
member that the input beam was unpolarized, so that the relative 7r/2 phase between field 
components was acquired by passing through the circular polarizer). Once again, the cor- 
rections are of the first order in the laboratory frame and second order in the beam frame. 
To conclude, note that about half of the intensity of the input beam was lost by passing 
across the polarizer: 



trJL = 1 ^trJ' 
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